Quantum lattice fluctuations in a frustrated Heisenberg spin-Peierls chain* 
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As a simple model for spin-Peierls systems we study a frustrated Heisenberg chain coupled to 
optical phonons. In view of the anorganic spin-Peierls compound CuGeOs we consider two different 
mechanisms of spin-phonon coupling. Combining variational concepts in the adiabatic regime and 
perturbation theory in the anti-adiabatic regime we derive effective spin Hamiltonians which cover 
the dynamical effect of phonons in an approximate way. Ground-state phase diagrams of these 
models are determined, and the effect of frustration is discussed. Comparing the properties of the 
ground state and of low-lying excitations with exact diagonalization data for the full quantum spin 
phonon models, good agreement is found especially in the anti-adiabatic regime. 

PACS numbers: 75.10.Jm, 63.20.Kr 
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I. INTRODUCTION 

The effect of a Peierls instability in quasi-one- 
dimensional spin systems, i.e. the instability of a uni- 
form spin chain towards dimerization induced by the in- 
teraction with lattice degrees of freedom, has attracted 
considerable attention over the last decades. Starting in 
the seventies with organic compounds of the TTF and 
TCNQ family [1], the interest in the Peierls instability 
was renewed with the discovery of a spin-Peierls (SP) 
transition in the anorganic compound CuGe03 in 1993 
by Hase et al. [2]. The most significant feature distin- 
guishing CuGeOs from other SP-compounds is the high 
energy of the involved optical phonons, which is compa- 
rable to the magnetic exchange integral J. In contrast to 
the organic materials no softening of these phonon modes 
is observed near the transition. Therefore the adiabatic 
treatment of the phonon subsystem used in the works of 
Pytte [3] or Cross and Fisher [4], does not seem appro- 
priate to describe the SP transition in CuGeO-3, although 
there are recent efforts in this direction [5]. Rather one 
has to take into account the effect of quantum lattice fluc- 
tuations which tend to decrease the SP transition tem- 
perature respectively the energy gap between the ground 
state and lowest excitations in the dimerized phase. Un- 
fortunately there are practically no analytic methods to 
handle coupled systems of spins (electrons) and phonons 
when all energy scales and coupling strengths are of the 
same order of magnitude. This is why many studies in- 
volving dynamical phonons rely on numerical methods, 
such as exact diagonalization (ED) [6,7], density matrix 
renormalization group (DMRG) [8,9] or Monte Carlo sim- 
ulation (MC) [10]. Only recently Zheng [11] developed an 
analytical approach to describe the SP instability of a XY 



spin chain, which is based on the unitary transformation 
method. It works well in the adiabatic and in the anti- 
adiabatic regime. In the latter case there are also some 
approaches to the Heisenberg spin chain interacting with 
optical phonons: Kuboki and Fukuyama [12] used per- 
turbation theory to derive an effective spin Hamiltonian, 
while Uhrig [13] applied the flow-equation method [14]. 

As a simple model which contains all important fea- 
tures of a SP system in the following we consider an 
antiferromagnetic Heisenberg chain coupled to a set of 
Einstein oscillators, 



with 



H — H s + H p + H sp 



H s = J 2j(S»-St+i + aS l -S l+2 ) , 

i 

Hp = cjq bfbj . 



(1) 

(2) 
(3) 



'Dedicated to H. Biittner on the occasion of his 60th 
birthday. 



The interaction of spins and phonons, H sp , can be mod- 
elled in two different ways, 

< c = ^(ft, + +&,)S,.S I+1 , (4) 

i 

H sp S =gJ2(K + bi)(Si-S i+1 - S r S^i) , (5) 

i 

where Si denote spin-i operators at lattice site i, while 
and bi are phonon creation and annihilation operators, 
respectively, Hl° c and differ in the mechanism, how 
the lattice influences the exchange integral. For Hl° c , the 
local coupling, one can think of a single harmonic degree 
of freedom directly modifying the magnetic interaction. 
In the context of CuGeC>3 this could correspond to side 
group effects (by the Germanium atoms) as discussed 
in Ref. [15]. In the case of , the difference coupling, 
the exchange depends directly on the spatial distance be- 
tween neighbouring spins. Note that it is not possible to 
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uniformly decrease or increase all exchange integrals with 
this type of spin phonon interaction. 

Although Hl° c seems to be more appropriate for 
CuGeC>3 we will consider both variants and compare 
their properties. In addition we take into account a frus- 
trating next-nearest neighbour interaction Ja, which in 
view of CuGeC>3 was introduced to explain susceptibility 
data [29] . As we will see below, the spin-phonon interac- 
tion is able to induce this kind of long ranged exchange 
as well. 

Motivated by the success of methods combining uni- 
tary transformations with variational and numerical 
techniques, which we used to study the Peierls transi- 
tion in the Holstein model of spinless fermions [16], and 
inspired by the work of Zheng [11], in this article we 
analyse the model (1) within the same framework. In 
particular we focus on the ground-state phase diagram 
as a function of spin-phonon coupling, frustration and 
phonon frequency, and compare our results with exact 
diagonalization data. 



II. EFFECTIVE SPIN MODELS 

To describe a static lattice dimerization in the adia- 
batic case of small phonon frequency ojq we start with a 
unitary transformation of H which shifts the equilibrium 
position of each oscillator by a constant amount alternat- 
ing from site to site, H = cxp(Si)H cxp(— Si) with 



si = |^£(-im + -M- 

For the terms involving phonons this yields 



(6) 



H n = 



+Nuj 



25 



H. 



loc 



diflf 



if>° c -A 7r £(-l) i S r S i+1 , 

i 

F s d f-2A^(-l) i S i -S i+1 . 



(7) 
(8) 
(9) 



| re) (n\H sp \m) (m\ 

En E m 



(11) 



where |re) and \m) denote eigenstates of the unperturbed 
Hamiltonian H s + H p with the corresponding eigenener- 
gies E n and E m . Since we do not know the exact eigen- 
states of H s , equation (11) can not be written in a simple 
form. One first guess which describes the phonon-part of 
52 in a correct way and ensures 5*2 to be antihermitean 
is 



s 2 oc =/-y>+-fo 4 )s 4 .s 4+1 , 

U)Q 



(12) 



S diff = f y_ yv 6 + _ 6 .)( Si .s i+1 - Si-Si_i) . (13) 

Wo 

Comparing with (11) loq gives the appropriate contribu- 
tion to the energy denominator, since _ff sp connects only 
phonon states differing by one in their occupation num- 
ber. What is missing is the contribution of spin states 
which have overlap via H sp . Therefore, using a numerical 
procedure we fix the free parameter / by the condition, 
that the amplitude of the state resulting from the appli- 
cation of H sp + [5*2, H s + H p 



pJ iA„=o to the ground state 
Figure 1 illustrates the variation 
of / with varying phonon frequency luq. We find that 



of H s + H p is minimal 
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Ajr will act as the variational parameter describing the 
dimerization of the system. 

Next we want to decouple spin and phonon degrees 
of freedom. The idea is to apply another unitary trans- 
formation H — exp(S , 2)-ff exp(— £2), where S2 should be 
chosen such that all contributions of first order in the 
coupling constant g disappear in the transformed Hamil- 
tonian. This gives the usual condition for a Schrieffer- 
Wolff transformation [17], 



H sp + [S2, H s + Hp] 



= 



(10) 



or 



FIG. 1.: Variation of / in the case of local coupling and 
lattice size N = 16. 

the general shape of / = f(u>o) depends only weakly on 
both, system size N and frustration a. While / — ► 1 in 
the anti-adiabatic frequency range (coq 3> J), the trans- 
formation exp(S l 2) vanishes completely as the frequency 
becomes small (u>o <C J). 

In contrast to electron-phonon systems with Holstein 
coupling, where a transformation similar to cxp(S , 2 ) com- 
pletely removes the electron-phonon interaction term, ap- 
plying the unitary transformation exp(S l 2) to H, we ob- 
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tain an infinite series of terms, which can not be summed 
up to a simple expression, i.e. 



(14) 



where [S 2 ,H]k denotes the iterated commutator 
[S 2 ,H] k+1 = [S 2 , [S 2 ,H} k ] with [S 2 ,H] = H. 

To get an effective spin model we consider only con- 
tributions up to fourth order in g and average over 
the phonon subsystem. The resulting spin Hamiltonian 
H c g = (H) contains long ranged Heisenberg interactions 
as well as numerous four- and six-spin couplings of the 
form (Sj-Sj)(Sfc-Sj) . . . (S m -S„). To a good approxima- 
tion we can neglect them and obtain 



(15) 



H cS /J = J + J2( J i + (-l)^)SrS i+ i 

i 

4 

i+k ■ 

i k=2 

Note that all phonon dynamics disappeared from the 
Hamiltonian £f e ffi but the effect of the spin-phonon in- 
teraction enters through both, the static dimcrization pa- 
rameter A^ and the different long range spin interactions. 

For the local coupling the corresponding interaction 
strengths are 



^loc _ _ 



Att 

J 



(16) 



J, 



loc 



N 



J 



1 

4A 



J 



j{ oc = 



8 Jy 2 ! 16 v 4 J 



rloc 
J 2 



rloc 

Jo — — 



+ -^"2 J - ( 4 ) + 
fg 2 Y(l - 2a) 



f 4 9 4 Y 2 
96 



(28 - 37a) , 



= a — 



•^Al -L)- / 4 9 4 F 2 -(1 - 2a) , 
4 v 4 y 96 v 



f 2 g 2 Ya , /V^ 



96 



(9 - 46a) , 



fr = 



9f 4 g 4 Y 2 a 
9(3 



while for the difference coupling we find 
S dis = 2 s ioc ^ 

LU 



rdiff 
J 



N 



J 

3\g 2 Y.f 
16 



(1 - I) 



(17) 
(18) 

(19) 

(20) 
(21) 

(22) 
(23) 



jr = l + 2A/(l-^) + 3/V 7-^ (24) 
+ g (1 - 4) + 24 ( 59 - 75a ) . 

12 v A' 24 v ; 
jf* = _/yy a _^£(l_Z) 

+ LJL ^ — (32- 119a), 



(25) 



(26) 



rdiff _ 

J A — 



48 

21/4 g 4y2 a 
48 



(27) 



To point out the relevant model parameters we intro- 
duced the dimensionless coupling constants A = g 2 / (Jojq) 
(cf. Refs. [3,4]) and g = g/ojo- Besides we used 



Y:=((b+-b l ) 2 ) = 



1 if T = 

coth(^) ifT>0 



(28) 



as a shorthand notation. To compare our result with that 
of Uhrig [13], in H d ff we have to set / = 1, which corre- 
sponds to the anti-adiabatic regime, and = 0. Indeed 
we recover (except for a prefactor 1/2 which in Ref. [13] 
enters erroneously going from eq. (11c) to eq. (13), and 
which consequently is also wrong in Ref. [9, eq. (3-4)]) 
all second order terms derived with the flow-equation 
method, supplemented by some new fourth order con- 
tributions. 

However, in contrast to Uhrig [13], who proposed to 
take a thermal average within the phonon subsystem, we 
believe, that averaging over the phonon vacuum, i.e. set- 
ting T — 0, is a more natural choice, especially if we 
are dealing with phonon-energies in the anti-adiabatic 
regime, where the interesting energy- and temperature- 
scales are small in respect of luq- Therefore in the follow- 
ing we use Y = — 1 and (bfbi) = exclusively. 



III. TRANSITION TO A GAPPED PHASE 

A prominent feature associated with the SP instability 
is of course the existence of an energy gap between the 
ground state and lowest excitations. Considering, in a 
first step, the pure spin model H s , it is known that the 
spectrum is gapless for the Heisenberg chain with a = 0, 
where the lowest spinon excitations (triplet and singlet) 
are degenerate with the ground state at momenta q = 
and 7r [18-20]. In contrast the system has a two-fold 
degenerate ground state and a gap to lowest triplet exci- 
tations at a = 0.5, the Majumdar- Ghosh point [21]. At 
some intermediate frustration a c the model undergoes a 
transition from the gapless to the gapped phase, which 
is of Kosterlitz-Thouless type [22-24] . Using arguments 
of conformal field theory one can show that the lowest 



3 



singlet and triplet excitations of a finite system of size N 
become degenerate at a c (N), where the dependence on 
N is only weak and a c (N) - a c (oo) ~ N~ 2 [25-27]. This 
was used in Refs. [28,29] to determine a c = 0.2411. 

Looking at our effective spin models H e s we find that 
the interaction with the optical phonons induces the same 
kind of frustrating next-nearest neighbour interaction. 
Therefore, without any explicit frustration a, the effec- 
tive frustration a c g := J2/J1 due to the phonons can 
lead to a gap in the energy spectrum and to spontaneous 
dimerization, as was already discussed in Ref. [12]. This 
effect is most important in the anti-adiabatic frequency 
range. 

Another mechanism producing a gap is (static) dimer- 
ization, i.e. an alternation of the nearest neighbour ex- 
change integral. Taking the adiabatic limit of our effec- 
tive model, / — ► and 5 — > A n , the ground-state energy 
of H s + 5^ i (-l) l Sj-Sj + i is known to deviate from its 
value at 5 = like <5 4 / 3 [4], while the elastic energy in- 
creases with 8 2 . Therefore for all couplings A the ground- 
state energy of H e g(f — 0) is minimal, if S is finite. At 
the same time proportional to (5 2 / 3 a gap opens in the 
spectrum. 

By taking into account both mechanisms we can now 
determine the transition from the gapless phase to the 
gapped one. As the SP-system behaves different for the 
two couplings, we treat them separately, starting with 
the local coupling case. 



A. Local coupling 

In a first step we set A w = and use the level-crossing 
criterium [28,29,9] to calculate the critical line in the a-g- 
plane for different phonon frequencies ojo and system sizes 
N. Since contains longer ranged interactions such 
as Si-Si + 3 and Si-Si+4, this line slightly deviates from 
the line a e s = ot c , and we have to calculate it separately. 
Applying the Lanczos algorithm to the effective model 
we obtain the critical line with high accuracy on local 
workstations (N < 20) . On the other hand we determine 
the level-crossing in the original model (1) by using the 
methods described in Ref. [6]. In the case wq/J = 0.1, the 
latter is complicated for the small spin phonon systems 
we can handle with our Lanczos diagonalization code on 
present day parallel computers, since the finite size gap 
to the first triplet excitation is a few times larger than wo- 
lf g is small enough, each eigenvalues of the spin system 
is accompanied by a ladder of phonon states. Therefore 
several 'copies' of the singlet ground state lie within the 
singlet-triplet gap if u>o -C J, and the singlet we consider 
for the level-crossing is not easy to determine unless we 
go to very large systems (see Ref. [9, Tab. II]). 

Figure 2(a) and 2(b) show the critical lines in the effec- 
tive (bold solid) and the original model (bold dashed) as 
well as the lines a e ff = a c (thin solid). As in the case of 
the pure spin model, the critical lines depend only weakly 
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FIG. 2.: Singlet-triplet level-crossing (solid lines) and on- 
set of dimerization (dashed with symbols) in the effective 
model in comparison to the level-crossing in the original 
model (bold dashed lines) at ujq/J = 10, 1 and 0.1. 



on N. We can therefore compare exact data for the orig- 
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inal model and N — 8 with data for the effective model 
and N = 16. While the results differ noticeable for in- 
termediate phonon frequency loq ~ 1, the agreement is 
excellent in the anti-adiabatic frequency range uo 3> J. 
With increasing wo the critical curve exhibits a remark- 
able upturn before crossing the abscissa, i.e. the frustra- 
tion is suppressed for small spin phonon coupling, but 
overcritical for strong coupling. It is this feature which 
makes it necessary to expand (14) up to fourth order to 
approximate H ioc in a correct way. A second order the- 
ory is not capable to describe the observed critical line. 

Another point we can study within our effective model 
is the behaviour of the critical spin phonon coupling g c 
at a = in the limit ujq/J — > co, i.e. the limit of the 
crossing point of the critical line and the abscissa. As the 
effects of the longer ranged interactions are rather small, 
we can solve the equation a l °g — a c . Setting / = 1 
(compare Fig. 1) and a = we find 

*=£ + ite) a+ f w,,h (29 > 

f^l-^lto.), (30) 

w 3,l . .37 7 . ,„„. 
Q = — -(- +a c ) - ( 1 a c ), (31) 

and in the limit of infinite phonon frequency g c ap- 
proaches a finite value, 

wo/J^oo y 3(1 + 2a c ) 
for the model with local coupling. 
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FIG. 3.: Low-lying excitations in the effective (filled sym- 
bols) and the original model (open symbols) for different 
frequency and local coupling g. 



In the case small phonon frequency, ujq -C J, the sec- 
ond transformation exp(S , 2 ) looses their importance, and 
the effective frustration due to the spin phonon interac- 
tion is replaced by the dimcrization as the relevant mech- 
anism leading to an energy gap. We account for this effect 
by allowing for a finite A w in our approximation. Using 
the Hcllmann-Feynman theorem and numerical diagonal- 
ization of finite spin systems we determine A v such that 
the ground-state energy of is minimal. Depend- 

ing on coupling strength g and frustration a the sys- 
tem prefers to remain in the undimerized, gapless phase 
(A ff = 0) or to develop a non-zero dimcrization leading 
to a gap. In Figs. 2(a) - (c) we plotted these transi- 
tion lines (dashed line with symbols) in addition to those 
obtained by level-crossing. As we already found in our 
study of the Holstein model of spinless fermions [16], for 
small luq the transition to the dimerized phase depends 
noticeable on the system size N (see Fig. 2(c)), while 
the finite-size dependence is weak in the anti-adiabatic 
regime (cf. Fig. 2(a)). In addition, for luq/J = 10 the 
transition is consistent with the critical line determined 
via level-crossing. 




FIG. 4.: Critical coupling g c versus frequency for if^r 
with a = 0. 

To compare properties of the original and the effec- 
tive model also in the case of small phonon frequency, 
we consider the dimcrization. As a quantity which cor- 
responds to <5 loc / diff we take the static (lattice) structure 
factor [6,7], 




j.k 



calculated in the ground state of (1), where Uj = b~j + 
bj. The inset of Fig. 2(c) displays {S/uj ) 2 for a = 0.36, 
luq/J = 0.1 (solid line), 0.316 (dashed line), and N = 8. 
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The results for the effective model (bold lines) and the 
original model (thin lines with symbols) agree rather well, 
especially for u>o/J = 0.1. 

Another feature we can compare is the dispersion of 
low-lying excitations. Figure 3 shows the energy of the 
lowest triplet excitations, calculated exactly and within 
our approximation. Clearly for uq/J = 0.3 the correct 
dispersion is flattened at momenta near q = it/2. This re- 
sults from the energy of the dispcrsionless phonons, which 
is added to the lowest triplet at q = n. Of course our ef- 
fective model does not contain these low-lying phonon 
excitations. However, as soon as luq > J the lowest ex- 
citations are due to renormalized spin interactions and 
well approximated by the effective model. 

To collect the results of this subsection we show in 
Fig. 4 the critical coupling g c (a — 0) over a wide range 
of phonon frequencies, using both criteria for the phase 
transition. Symbols stand for the onset of dimerization, 
while the bold line corresponds to a c g = a c . As ex- 
pected, we find that our approximation is somewhat un- 
reliable for intermediate phonon frequency. The singular- 
ity of g c at uj / J — 1 is a manifestation of this deficiency. 
The correct critical line will connect adiabatic and anti- 
adiabatic behaviour in a continuous way (compare also 
next subsection and [9]). 



B. Difference coupling 

The procedure to determine the phase transition in 
the SP-system with difference coupling is the same as 
described before. In the anti-adiabatic regime we set 
= and calculate the position of the crossing of 
the first triplet and the first singlet excitation for both, 
the original and the effective model. The results for 
u>o/J = 10 and 1 arc shown in Figs. 5(a) and 5(b) respec- 
tively. In contrast to the local coupling the structure of 
the critical line for high phonon frequency (u>o/J = 10) 
is much simpler. It appears that one would get the same 
shape also for a second order theory. However, to enlarge 
the application area of our approximation taking into ac- 
count higher order contributions is still appropriate. As 
before, the agreement between the original and the ef- 
fective model is excellent in the anti-adiabatic regime, 
while the deviations increase with approaching interme- 
diate frequencies. 

Calculating the behaviour of the critical coupling, 
g c {ct = 0) in the limit of infinite phonon frequency, 
loq/ J — > oo, we now find 
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FIG. 5.: Level-crossing (solid lines) and onset of dimer- 
ization (dashed with symbols) in the effective model in 
comparison to the level-crossing in the original model 
(bold dashed lines) at ujq/J = 10, 1 and 0.1. 



and, different to the local coupling case, g c tends to zero, 



lim g c = . 

wo/ J— >oo 



(37) 



G 



While the q = and the q — ir phonon mode compete in 
the case of local spin phonon coupling, allowing for a sta- 
ble gapless phase up to a critical g, there is no interaction 
with the q = mode in . Therefore the q = tt mode 
induces long ranged exchange more efficiently, leading to 
a vanishing g c for u; /J^> oo. 

For small phonon frequencies wo -C J again we deter- 
mine the optimal dimerization A T and the critical line 
beyond which A w starts to be non-zero. The results 
are shown in Figs. 5(b) and 5(c). For large phonon fre- 
quency we find an unstable behaviour of A n , it is finite 
for some <?, but vanishes before growing to substantial 
values again. We therefore make no attempt to fix the 
onset of ^ for cj / J > 1 . 

As in the case of local coupling the dimerization in the 
original and the effective model agrees well for uj /J = 
0.1 (inset Fig. 5(c)). 
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FIG. 6.: Critical coupling g c versus frequency for 
with a = 0. 

Finally in Fig. 6 we combine the above results to ob- 
tain the phase diagram in the J/uoQ-g plane. Except for 
wo/ J — ► oo, where g c — ► 0, the behaviour is similar to the 
case of local coupling. Again at loq/J ~ 1 the effective 
model gives a singularity in g c , while the correct result 
should be continuous. Comparing the transition line with 
the very recent DMRG data of Bursill ct al. [9], we find 
very good agreement in the anti-adiabatic regime (the 
phase transition does not change, going from the second 
to the fourth order effective theory). For small phonon 
frequency loq/ J — > the results of Bursill et al. suggest a 
finite limit for g c /u>o, while our model gives an increasing 
ratio of g c /^o- However, as the finite size effects are large 
for luq <C J, determining the correct value of g c is a del- 
icate procedure and the exact g c might be much smaller 
than depicted in Fig. 6. 



IV. GROUND-STATE PHONON DISTRIBUTION 

In the course of exact diagonalizations of the phonon 
dynamical model (1) we observed another interest- 
ing feature distinguishing the two mechanisms of spin 
phonon interaction. Turning our attention to the 
phonon distribution in the ground state of (1), we find 
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FIG. 7.: Even/odd-phonon distribution in the ground 
state of the model with difference coupling (g = 1.4, 
ojq/J = 0.1) and mechanisms that suppress the effect: 
local coupling, frustration and anti-adiabaticity. 

for the model with difference coupling in the adia- 
batic regime, that the system prefers states with even 
phonon occupation numbers (cf. Fig. 7(a)). This be- 
haviour reminds of a simple two-level system, more pre- 
cisely the Rabi (pseudo-Jahn- Teller) Hamiltonian [30] 
H = Aa z + g(b + + b)a x + ujob + b <7o, where an interac- 
tion with a bosonic degree of freedom connects the lev- 
els. Solving this model in the adiabatic strong-coupling 
case, uiq/ A -C 1 and g = g/uJo > 1, we obtain a similar 
even/odd oscillation in the ground-state phonon distri- 
bution. One can understand this effect within standard 
perturbation theory. Starting with the two levels ±A and 
the corresponding eigenstates (J) and (J) , it is obvious, 
that adding only one phonon requires an extra amount 
of 2 A exchange energy, because the system is excited to 
the upper level. Therefore the system prefers an even 
number of phonons for the ground-state wave function. 
In the anti-adiabatic case, ujq » A, to a good approxi- 
mation we can omit the Aa z term and solve the model by 
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unitary transformation, leading to a Poisson distribution 
of the phonons (coherent state). 

In view of the SP system with difference coupling, the 
q = 7r phonon is the most relevant (in fact in the ground 
state almost all phonons occupy this mode), and adding 
one phonon changes the momentum of the spin system 
by 7r. Consequently the gap between the (singlet) ground 
state at q = (N = 4k, k e N) or ir (N = 4k + 2) and the 
lowest singlet state at q = it, or respectively, plays the 
role of 2A in the two-level system. As this singlet-singlet 
gap is only due to the finite size of the considered spin 
chain, we believe that the even/odd oscillations will dis- 
appear in the thermodynamic limit. The reduction (dis- 
appearance) of the even/odd imbalance resulting from a 
finite frustration with a much smaller singlet-singlet gap 
(zero at a = 0.5), cf. Fig. 7(c), can be taken as a first 
indication. Of course we observe a smooth phonon dis- 
tribution in the anti-adiabatic case (Fig. 7(d)). 

The model with local spin phonon coupling exhibits 
the usual Poisson distribution for all frequencies, since 
the interaction is due to the q — it and the q = phonon 
mode (Fig. 7(b)). 



phase transition by means of the static dimerization A n , 
which changes from zero at small spin phonon coupling 
to a finite value beyond a critical coupling. 

At intermediate frequency the situation remains unsat- 
isfactory as the critical coupling behaves discontinuously. 
Here numerical methods, including the full phonon dy- 
namics, still provide the only reliable tool to study the 
transition. Nevertheless the proposed effective models 
help to understand the physical mechanisms leading to 
spontaneous dimerization of the interacting spin-phonon 
system. 
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V. CONCLUSION 



In summary, we have studied the spin-Pcicrls insta- 
bility of a frustrated Heisenberg spin chain coupled to 
optical phonons of energy o>o- Using the concept of uni- 
tary transformations we derive effective spin Hamilto- 
nians, which cover the spin phonon interaction by two 
mechanisms, static dimerization A T and long ranged ex- 
change couplings. Both can lead to an energy gap be- 
tween the ground state and lowest excitations, which is 
related to a Peierls instability of the spin system. In 
the anti-adiabatic phonon frequency range, uio ^> J, we 
verify and extend the Hamiltonian obtained recently with 
the flow-equation method [13], while we recover the usual 
static SP model in the limit loq/J — > 0. 

To determine the transition to the gapped phase in 
the case of large phonon frequency we use the level- 
crossing criterium, which proved to be very accurate 
for similar models [28,29,9]. For the two types of spin 
phonon coupling {local: U{ Sj-Sj+i, and difference: {ui — 
Ui + i) Sj-Sj+i) we consider here, the results of our effec- 
tive models agree very well with data from exact diago- 
nalization of the original, phonon dynamical model and 
with recent DMRG data [9] (difference coupling only). 
In the case of local coupling two phonon modes (q = 
and 7r) compete and allow for a gapless phase to exist 
in a wider parameter range. Furthermore we observe a 
non-monotonous behaviour of the phase transition line as 
a function of spin-phonon interaction and frustration a. 
With difference coupling spins and phonons interact al- 
most only through the 7r-mode, which is able to dimerize 
the system most efficiently. 

For phonon frequencies ujq < J, we determine the 
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